interpolation in numerical method